Patterns of Recent Natural Selection on Genetic Loci Associated with Sexually Differentiated Human Body Size and Shape Phenotypes

Levels of sex differences for human body size and shape phenotypes are hypothesized to have adaptively reduced following the agricultural transition as part of an evolutionary response to relatively more equal divisions of labor and new technology adoption. In this study, we tested this hypothesis by studying genetic variants associated with five sexually differentiated human phenotypes: height, body mass, hip circumference, body fat percentage, and waist circumference. We first analyzed genome-wide association (GWAS) results for UK Biobank individuals (~194,000 females and ~167,000 males) to identify a total of 114,199 single nucleotide polymorphisms (SNPs) significantly associated with at least one of the studied phenotypes in females, males, or both sexes (P \u3c 5x10-8). From these loci we then identified 3,016 SNPs (2.6%) with significant differences in the strength of association between the female- and male-specific GWAS results at a low false-discovery rate (FDR \u3c 0.001). Genes with known roles in sexual differentiation are significantly enriched for co-localization with one or more of these SNPs versus SNPs associated with the phenotypes generally but not with sex differences (2.73-fold enrichment; permutation test; P = 0.0041). We also confirmed that the identified variants are disproportionately associated with greater phenotype effect sizes in the sex with the stronger association value. We then used the singleton density score statistic, which quantifies recent (within the last ~3,000 years; post-agriculture adoption in Britain) changes in the frequencies of alleles underlying polygenic traits, to identify a signature of recent positive selection on alleles associated with greater body fat percentage in females (permutation test; P = 0.0038; FDR = 0.0380), directionally opposite to that predicted by the sex differences reduction hypothesis. Otherwise, we found no evidence of positive selection for sex difference-associated alleles for any other trait. Overall, our results challenge the longstanding hypothesis that sex differences adaptively decreased following subsistence transitions from hunting and gathering to agriculture

The Kabua 1 cranium: Virtual anatomical reconstructions

Our current knowledge of the emergence of anatomically modern humans, and the human lineage in general, is limited, in large part because of the lack of a well preserved and well dated fossil record from Pleistocene Africa. Thus, the primary aim of our research is to partly relieve this problem by virtually reconstructing and analyzing the hominin cranial remains of Kabua 1, found in Kenya in the 1950s. Most scholars have argued that Kabua 1 represents an anatomically modern Homo sapiens, although the fragmentary nature of the remains and lack of a chronometric date hinder robust phylogenetic and taxonomic assessments. This manuscript presents the first steps taken to resolve this issue, namely a set of reconstructions of the specimen that would allow comparison with the fossil record. First, we virtually removed sediment and laboratory adhesives from μct scans of the fragments. Subsequently, all fragments were separated by segmentation of the μct data and described. Finally, virtual surface projections were used in the creation of several anatomical reconstructions, based on separate reference crania. These first steps provide a framework that will be used for quantitative shape analyses that aim to more firmly place these remains in the context of human evolution

Patterns of recent natural selection on genetic loci associated with sexually differentiated human body size and shape phenotypes.

Levels of sex differences for human body size and shape phenotypes are hypothesized to have adaptively reduced following the agricultural transition as part of an evolutionary response to relatively more equal divisions of labor and new technology adoption. In this study, we tested this hypothesis by studying genetic variants associated with five sexually differentiated human phenotypes: height, body mass, hip circumference, body fat percentage, and waist circumference. We first analyzed genome-wide association (GWAS) results for UK Biobank individuals (~194,000 females and ~167,000 males) to identify a total of 114,199 single nucleotide polymorphisms (SNPs) significantly associated with at least one of the studied phenotypes in females, males, or both sexes (P<5x10-8). From these loci we then identified 3,016 SNPs (2.6%) with significant differences in the strength of association between the female- and male-specific GWAS results at a low false-discovery rate (FDR<0.001). Genes with known roles in sexual differentiation are significantly enriched for co-localization with one or more of these SNPs versus SNPs associated with the phenotypes generally but not with sex differences (2.73-fold enrichment; permutation test; P = 0.0041). We also confirmed that the identified variants are disproportionately associated with greater phenotype effect sizes in the sex with the stronger association value. We then used the singleton density score statistic, which quantifies recent (within the last ~3,000 years; post-agriculture adoption in Britain) changes in the frequencies of alleles underlying polygenic traits, to identify a signature of recent positive selection on alleles associated with greater body fat percentage in females (permutation test; P = 0.0038; FDR = 0.0380), directionally opposite to that predicted by the sex differences reduction hypothesis. Otherwise, we found no evidence of positive selection for sex difference-associated alleles for any other trait. Overall, our results challenge the longstanding hypothesis that sex differences adaptively decreased following subsistence transitions from hunting and gathering to agriculture

A genome-wide association scan implicates <i>DCHS2, RUNX2, GLI3, PAX1</i> and <i>EDAR</i> in human facial variation

We report a genome-wide association scan for facial features in ∼6,000 Latin Americans. We evaluated 14 traits on an ordinal scale and found significant association (P values−8) at single-nucleotide polymorphisms (SNPs) in four genomic regions for three nose-related traits: columella inclination (4q31), nose bridge breadth (6p21) and nose wing breadth (7p13 and 20p11). In a subsample of ∼3,000 individuals we obtained quantitative traits related to 9 of the ordinal phenotypes and, also, a measure of nasion position. Quantitative analyses confirmed the ordinal-based associations, identified SNPs in 2q12 associated to chin protrusion, and replicated the reported association of nasion position with SNPs in PAX3. Strongest association in 2q12, 4q31, 6p21 and 7p13 was observed for SNPs in the EDAR, DCHS2, RUNX2 and GLI3 genes, respectively. Associated SNPs in 20p11 extend to PAX1. Consistent with the effect of EDAR on chin protrusion, we documented alterations of mandible length in mice with modified Edar function

Ancient Jomon genome sequence analysis sheds light on migration patterns of early East Asian populations

Funder: The excavation of the Ikawazu Jomon individual was supported by Grant-in-Aid for Scientific Research (B) (25284157) to YY. The Ikawazu Jomon genome project was organized by HI, and TH & HO who were supported by MEXT KAKENHI Grant Numbers 16H06408 and 17H05132, by Grant-in-Aid for Scientific Research on Innovative Areas (Cultural History of Paleoasia), and by Grant-in-Aid for Challenging Exploratory Research (23657167) and for Scientific Research (B) (17H03738). The Ikawazu Jomon genome sequencing was supported by JSPS KAKENHI Grant Number 16H06279 to ATo, and partly supported in the CHOZEN project in Kanazawa University, and in the Cooperative Research Project Program of the Medical Institute of Bioregulation, Kyushu University. Computations for the Ikawazu Jomon genome were partially performed on the NIG supercomputer at ROIS National Institute of Genetics.Abstract: Anatomically modern humans reached East Asia more than 40,000 years ago. However, key questions still remain unanswered with regard to the route(s) and the number of wave(s) in the dispersal into East Eurasia. Ancient genomes at the edge of the region may elucidate a more detailed picture of the peopling of East Eurasia. Here, we analyze the whole-genome sequence of a 2,500-year-old individual (IK002) from the main-island of Japan that is characterized with a typical Jomon culture. The phylogenetic analyses support multiple waves of migration, with IK002 forming a basal lineage to the East and Northeast Asian genomes examined, likely representing some of the earliest-wave migrants who went north from Southeast Asia to East Asia. Furthermore, IK002 shows strong genetic affinity with the indigenous Taiwan aborigines, which may support a coastal route of the Jomon-ancestry migration. This study highlights the power of ancient genomics to provide new insights into the complex history of human migration into East Eurasia

Studies of new Higgs boson interactions through nonresonant HH production in the b¯bγγ fnal state in pp collisions at √s = 13 TeV with the ATLAS detector

A search for nonresonant Higgs boson pair production in the b ¯bγγ fnal state is performed using 140 fb−1 of proton-proton collisions at a centre-of-mass energy of 13 TeV recorded by the ATLAS detector at the CERN Large Hadron Collider. This analysis supersedes and expands upon the previous nonresonant ATLAS results in this fnal state based on the same data sample. The analysis strategy is optimised to probe anomalous values not only of the Higgs (H) boson self-coupling modifer κλ but also of the quartic HHV V (V = W, Z) coupling modifer κ2V . No signifcant excess above the expected background from Standard Model processes is observed. An observed upper limit µHH &lt; 4.0 is set at 95% confdence level on the Higgs boson pair production cross-section normalised to its Standard Model prediction. The 95% confdence intervals for the coupling modifers are −1.4 &lt; κλ &lt; 6.9 and −0.5 &lt; κ2V &lt; 2.7, assuming all other Higgs boson couplings except the one under study are fxed to the Standard Model predictions. The results are interpreted in the Standard Model efective feld theory and Higgs efective feld theory frameworks in terms of constraints on the couplings of anomalous Higgs boson (self-)interactions

Measurement of the H → γ γ and H → ZZ∗ → 4 cross-sections in pp collisions at √s = 13.6 TeV with the ATLAS detector

The inclusive Higgs boson production cross section is measured in the di-photon and the Z Z∗ → 4 decay channels using 31.4 and 29.0 fb−1 of pp collision data respectively, collected with the ATLAS detector at a centre of-mass energy of √s = 13.6 TeV. To reduce the model dependence, the measurement in each channel is restricted to a particle-level phase space that closely matches the chan nel’s detector-level kinematic selection, and it is corrected for detector effects. These measured fiducial cross-sections are σfid,γ γ = 76+14 −13 fb, and σfid,4 = 2.80 ± 0.74 fb, in agreement with the corresponding Standard Model predic tions of 67.6±3.7 fb and 3.67±0.19 fb. Assuming Standard Model acceptances and branching fractions for the two chan nels, the fiducial measurements are extrapolated to the full phase space yielding total cross-sections of σ (pp → H) = 67+12 −11 pb and 46±12 pb at 13.6 TeV from the di-photon and Z Z∗ → 4 measurements respectively. The two measure ments are combined into a total cross-section measurement of σ (pp → H) = 58.2±8.7 pb, to be compared with the Stan dard Model prediction of σ (pp → H)SM = 59.9 ± 2.6 p

Genomic analyses inform on migration events during the peopling of Eurasia.

High-coverage whole-genome sequence studies have so far focused on a limited number of geographically restricted populations, or been targeted at specific diseases, such as cancer. Nevertheless, the availability of high-resolution genomic data has led to the development of new methodologies for inferring population history and refuelled the debate on the mutation rate in humans. Here we present the Estonian Biocentre Human Genome Diversity Panel (EGDP), a dataset of 483 high-coverage human genomes from 148 populations worldwide, including 379 new genomes from 125 populations, which we group into diversity and selection sets. We analyse this dataset to refine estimates of continent-wide patterns of heterozygosity, long- and short-distance gene flow, archaic admixture, and changes in effective population size through time as well as for signals of positive or balancing selection. We find a genetic signature in present-day Papuans that suggests that at least 2% of their genome originates from an early and largely extinct expansion of anatomically modern humans (AMHs) out of Africa. Together with evidence from the western Asian fossil record, and admixture between AMHs and Neanderthals predating the main Eurasian expansion, our results contribute to the mounting evidence for the presence of AMHs out of Africa earlier than 75,000 years ago.Support was provided by: Estonian Research Infrastructure Roadmap grant no 3.2.0304.11-0312; Australian Research Council Discovery grants (DP110102635 and DP140101405) (D.M.L., M.W. and E.W.); Danish National Research Foundation; the Lundbeck Foundation and KU2016 (E.W.); ERC Starting Investigator grant (FP7 - 261213) (T.K.); Estonian Research Council grant PUT766 (G.C. and M.K.); EU European Regional Development Fund through the Centre of Excellence in Genomics to Estonian Biocentre (R.V.; M.Me. and A.Me.), and Centre of Excellence for Genomics and Translational Medicine Project No. 2014-2020.4.01.15-0012 to EGC of UT (A.Me.) and EBC (M.Me.); Estonian Institutional Research grant IUT24-1 (L.S., M.J., A.K., B.Y., K.T., C.B.M., Le.S., H.Sa., S.L., D.M.B., E.M., R.V., G.H., M.K., G.C., T.K. and M.Me.) and IUT20-60 (A.Me.); French Ministry of Foreign and European Affairs and French ANR grant number ANR-14-CE31-0013-01 (F.-X.R.); Gates Cambridge Trust Funding (E.J.); ICG SB RAS (No. VI.58.1.1) (D.V.L.); Leverhulme Programme grant no. RP2011-R-045 (A.B.M., P.G. and M.G.T.); Ministry of Education and Science of Russia; Project 6.656.2014/K (S.A.F.); NEFREX grant funded by the European Union (People Marie Curie Actions; International Research Staff Exchange Scheme; call FP7-PEOPLE-2012-IRSES-number 318979) (M.Me., G.H. and M.K.); NIH grants 5DP1ES022577 05, 1R01DK104339-01, and 1R01GM113657-01 (S.Tis.); Russian Foundation for Basic Research (grant N 14-06-00180a) (M.G.); Russian Foundation for Basic Research; grant 16-04-00890 (O.B. and E.B); Russian Science Foundation grant 14-14-00827 (O.B.); The Russian Foundation for Basic Research (14-04-00725-a), The Russian Humanitarian Scientific Foundation (13-11-02014) and the Program of the Basic Research of the RAS Presidium “Biological diversity” (E.K.K.); Wellcome Trust and Royal Society grant WT104125AIA & the Bristol Advanced Computing Research Centre (http://www.bris.ac.uk/acrc/) (D.J.L.); Wellcome Trust grant 098051 (Q.A.; C.T.-S. and Y.X.); Wellcome Trust Senior Research Fellowship grant 100719/Z/12/Z (M.G.T.); Young Explorers Grant from the National Geographic Society (8900-11) (C.A.E.); ERC Consolidator Grant 647787 ‘LocalAdaptatio’ (A.Ma.); Program of the RAS Presidium “Basic research for the development of the Russian Arctic” (B.M.); Russian Foundation for Basic Research grant 16-06-00303 (E.B.); a Rutherford Fellowship (RDF-10-MAU-001) from the Royal Society of New Zealand (M.P.C.)

Search for flavour-changing neutral tqH interactions with H → γγ in pp collisions at √s = 13 TeV using the ATLAS detector

A search for flavour-changing neutral interactions involving the top quark, the Higgs boson and an up-type quark q (q = c, u) is presented. The proton-proton collision data set used, with an integrated luminosity of 139 fb−1, was collected at √s = 13 TeV by the ATLAS experiment at the Large Hadron Collider. Both the decay process t → qH in tt¯ production and the production process pp → tH, with the Higgs boson decaying into two photons, are investigated. No significant excess is observed and upper limits are set on the t → cH and the t → uH branching ratios of 4.3 × 10−4 and 3.8 × 10−4, respectively, at the 95% confidence level, while the expected limits in the absence of signal are 4.7 × 10−4 and 3.9 × 10−4. Combining this search with ATLAS searches in the H → τ+τ− and H → bb¯ final states yields observed (expected) upper limits on the t → cH branching ratio of 5.8 × 10−4 (3.0 × 10−4) at the 95% confidence level. The corresponding observed (expected) upper limit on the t → uH branching ratio is 4.0 × 10−4 (2.4 × 10−4)